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Abstract 

We use the doped Pullerenes as an example of how reahstic systems can be described 
by simple models. Starting from the band structure we set up a tight-binding model 
that describes the ti^ conduction band. Adding correlation terms we arrive at a 
generalized Hubbard Hamiltonian that we treat using quantum Monte Carlo. To 
address the problem of superconductivity in the doped Fullerenes, we study the 
screening of a point charge. We find surprisingly efficient metallic screening even for 
strong correlations, almost up to the Mott transition, and discuss the implications on 
superconductivity, in particular the effect of the efficient screening on the Coulomb 
pseudopotential and the electron-phonon coupling. 
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1 Introduction 

The sheer size of the many-body Hilbert space makes treating strongly cor- 
related systems adequately extremely difficult or even impossible. Examples 
of considerable interest are the superconducting doped Fullerenes. Even for 
a single Cgo molecule a full many-body calculation is still a challenge and 
calculations for solids made of Fullerenes are simply out of question. In this 
situation we are forced to restrict our attention to only the most relevant 
degrees of freedom. For the doped Fullerenes these are the electrons in the 
ti„-band. Starting from ah initio density functional calculations we set up a 
tight-binding Hamiltonian that describes the electrons in the ti„-band only. 
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Including correlation effects we arrive at a generalized Hubbard Hamiltonian 
that can be treated hj quantum Monte Carlo (QMC). Using the fixed-node 
approximation we study the screening of an point charge by the electrons 
in the ti„-band. We find that the screening is surprisingly efficient even for 
strong correlations, almost up to the Mott transition. This has important im- 
plications for superconductivity in the doped FuUerenes. Given that molecular 
vibration energies are of the same order as electronic energies, retardation ef- 
fects are inefficient at reducing the electron-electron repulsion. It is therefore 
not clear how the weak electron-phonon attraction can lead to superconduc- 
tivity. Efficient metallic screening, as found in our calculations, can, however, 
reduce the electron-electron repulsion enough to allow for an electron-phonon 
driven superconductivity. But the screening does, of course, also affect the 
coupling to the phonons. It turns out that, due to screening, the alkali and Ag 
modes couple only weakly, while the couphng to the Hg modes is not affected. 
Therefore, although being electron-phonon driven, superconductivity in the 
doped Fullerenes differs in important ways from the conventional picture of 
superconductivity. 



2 Model Hamiltonian 

Treating correlations in the Fullerenes is quite difficult. Even for a single Cgo 
molecule a full many-body calculation is still a challenge and simulations of 
Fullerenes, which are solids made of Cgo molecules are well beyond current 
computational capabilities. Solid Ceo is characterized by weak inter-molecular 
interactions. Hence the molecular levels merely broaden into narrow, well sep- 
arated bands [1] . Doping the solid with alkali metals has the effect of filling the 
band originating from the molecular ti„-level with the weakly bound valence 
electrons of the alkali atoms. Since the ti„-level is three-fold degenerate the 
corresponding band can hold up to six electrons per molecule; for AsCeo, e.g., 
the band is half-filled. When we are interested in the low-energy properties 
of the tiu electrons, it is a good approximation to focus only on the region 
around the Fermi level; i.e. on the ti„-band, projecting out all the other bands 
[2]. That way we arrive at a tight-binding Hamiltonian comprising only the 
tiu orbitals, that reproduces the ab initio band structure remarkably well [3]. 

To obtain a realistic description we have to include the correlation arising 
from the Coulomb repulsion among the electrons on the same molecule. The 
resulting Hamiltonian describes the interplay between the hopping of electrons 
between different molecules and the Coulomb repulsion among the electrons 
on the same molecule 

(ij) nn'a i (na)<{n' a') 
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Here cj^^. creates an electron of spin a in the orbital with index n e {1, . . . , 3} 
on molecule i, the ti„j„' are hopping matrix elements between orbitals n and 
n' on neighboring molecules, and riina- = cl^crCina are occupation operators. 
Varying the value of the interaction term U, we can study the effect of cor- 
relations. The physical value for the doped FuUerenes is U ^ 1.2 — 1.4 eV [4], 
which has to be compared to the width of the ii„-band, W Ri 0.5 — 0.85 eV. 



3 Quantum Monte Carlo 

Wc now give a very brief outline of the quantum Monte Carlo method for 
determining ground states. The basic idea is to use the Hamiltonian to project 
out the ground-state from some trial function |^t), that we have guessed: 

e-^^|*T)'-=^|*o)- (2) 

To see how this works let us assume we knew the expansion of the trial function 
in eigenfunctions of H 

|*T) = Ecn|*n) e-^^|*T)=Ecne-^^M*n). (3) 

n 

Thus the component with the lowest eigenenergy is damped least during the 
projection, i.e., if cq 7^ 0, in the limit of large r the ground-state component 
will dominate. In practice we use for Hamiltonians H with a spectrum that is 
bounded, both, from below and above, an iterative projection of the form [5] 

1^(0)) = 1^^) |^(")) = [1 -r(if-«;)]|*('^-i)), (4) 

which, for small but finite r, also exactly projects to the ground-state. We see 
that the basic operation in (4) is a matrix- vector product. Since we are working 
in a many-body Hilbert space, the dimension of the vectors is, however, in 
general enormous; see Table 1 for an illustration. 

To understand the Monte Carlo method for doing the iteration we first rewrite 
(4) in configuration space. Here R denotes a configuration of the electrons in 
real space 

El^')(^'l^^"^) = E (i^ll - r{H - Eo)|i^K i^|^("-^)). (5) 



We see that the matrix F{R', K) maps configuration R into configurations 
R'. We clearly cannot follow every possible new configuration R' since that 
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Table 1 

Dimension of the Hilbert space for a system with + electrons on M lattice 
sites. For our model of doped Fullcrcncs AI = SA'^moi- The last column gives the 
amount of memory needed to store a single state vector. 

would lead to an exponential growth in the number of configurations as we 
iterate. The idea of Monte Carlo is then to sample only one of the configu- 
rations R' with a probability p{R\R). To do that we want to interpret the 
matrix elements of F{R',R) as probabihties. They are, however, in general 
not normalized and can even be negative 

F{R', R) = p{R', R) m{R', R) . (6) 
^ V ' ^ V ' 

prob. norm.&sign 

Normalization introduces the need for population control, while negative ma- 
trix elements introduce the sign problem [6,7]. 



4 Screening of a Point Charge 

In conventional superconductors the electron-phonon interaction leads to an 
effective electron-electron attraction. This attraction is, of course, counter- 
acted by the Coulomb repulsion between the electrons. In the conventional 
picture this repulsion is, however, strongly reduced by retardation effects [8]. 
The resulting effective Coulomb repulsion is described by the dimensionless 
Coulomb pseudopotential /i*, which is believed to be of the order of 0.1. For 
the doped Fullcrcncs the situation is different: Retardation effects are ineffi- 
cient [9]. Therefore the screening of the Coulomb interaction becomes impor- 
tant for reducing the electron-electron repulsion. Assuming that the random 
phase approximation (RPA) is valid for the electrons within the ti^-band it 
was found that efficient metallic screening significantly reduces the Coulomb 
pseudopotential [10]. In this scenario the Coulomb pseudopotential yU* ^ 0.3 
is substantially larger than that for conventional superconductors, but it is 
not too large to prevent superconductivity from being driven by the electron- 
phonon interaction. For strongly correlated systems like the doped PuUerenes 



4 



the use of the RPA seems, however, highly questionable. 

To address this question, we have investigated how well the RPA describes the 
screening within the fi„-band. It is clear that the RPA properly describes the 
screening when the kinetic energy is much larger than the interaction energy, 
i.e. RPA works well in the weakly correlated limit. For strong correlations, 
where the Coulomb energy dominates, the RPA gives qualitatively wrong re- 
sults: Introducing a test charge g on a molecule, RPA predicts that almost the 
same amount of electronic charge moves away form that molecule, since for a 
Coulomb integral U much larger than the band width W the gain in potential 
energy dominates the cost in kinetic energy. The RPA neglects, however, that 
in this limit, when an electron leaves a molecule it has to find another molecule 
with a missing electron, or there will be a large increase in Coulomb energy. 
It is not clear for what value of the Coulomb interaction U this qualitative 
breakdown of the RPA starts, and up to which values of U RPA gives still 
accurate results. To find out we have performed quantum Monte Carlo calcu- 
lations. To study the screening of a point charge q on molecule with index c 
we consider the Hamiltonian 

^ ~ $Z ^ tin,jn' C^inc^jn'cr (7) 
(ij) nn'a i (na)<(n'a') 

which differs from (1) only by the last term describing the interaction with 
the test charge. As a trial function we use a generalized Gutzwiller function 
of the form 

\^^)^g^g-^\^,), (8) 

where |$o) is a Slater determinant, D is the number of double occupancies 
in the system, and ric is the number of electrons on the molecule with the 
test charge. Calculating the expectation value nciVMC) = (\['t|^^c|^t) by 
variational Monte Carlo and the mixed estimator nc{DMC) = (\t'T|^c|^o) 
by fixed-node diffusion Monte Carlo, we obtain the ground-state expectation 
value Tie — (^ol^d^o) from the extrapolated estimator Uc ~ 2nc{DMC) — 
UciyMC). To estimate the accuracy of our approach we have compared the 
results of the quantum Monte Carlo (QMC) calculations with the results from 
exact diagonalization for a cluster of 4 molecules. We find that the QMC 
calculations are accurate up to very large values of the Coulomb interaction 
U. Performing QMC calculations for clusters of sizes A^moi — 32, 48, 64, 72, 
and 108 molecules, where exact diagonalization is not possible (cf. Table 1), we 
have extrapolated the screening charge Anc = 'n.c(O) — ndq) to infinite cluster- 
size, assuming a finite-size scaling of the form Aric = Anc(A^moi) + ct/A^moi 
[1 1] . The finite-size extrapolation gives only a small correction to the screening 
charge found for the larger clusters. Results for q — 0.25 e are shown in Figure 
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Fig. 1. Screening cliarge An on tlic site of a test charge g = 1/4 electron charges 
as a function of U /W, extrapolated to infinite cluster size. 

1. For small values of U the RPA somewhat underestimates the screening, 
a behavior also found in the electron gas [12]. For intermediate values of U 
{U/W ~ 1.0 — 2.0) the RPA still gives surprisingly accurate results, while for 
larger U it rapidly becomes qualitatively wrong, as discussed above. We thus 
find efficient, RPA-hke screening even for quite strong correlations close to the 
Mott transition. 



5 Screening and Electron-Phonon Coupling 

As pointed out in the preceeding discussion, efficient metallic screening helps 
to reduce the effective electron-electron repulsion, i.e. the Coulomb pseudopo- 
tential. But the screening also affects the electron-phonon coupling. At first it 
might appear that efficient screening is not really helpful for superconductiv- 
ity. Phonons couple to the electrons by perturbing the potential seen by the 
electrons. An example is the longitudinal mode of a jellium. Efficient screening 
tends to weaken the coupling to such phonons, since it reduces the perturba- 
tion considerably: The bare coupling constant g is reduced to g/e, where e is 
the dielectric constant [8]. To some extent, such a reduction is also at work 
in the FuUerenes. An example are the alkali phonons. Each Geo molecule is 
surrounded by 14 alkali ions that are bound with quite weak force constants. 
They should therefore respond strongly when an electron arrives on a Ceo 
molecule. This was, however, not confirmed experimentally; an alkali isotope 
effect could, e.g., not be observed [13]. Given the efficient metallic screening 
this finding can be naturally understood. When an electron arrives on a Geo 
molecule, other electrons leave the molecule, which thus stays almost neutral. 
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Fig. 2. Schematic picture of the change in the energy levels as function of the 
phonon coordinate Q (deformation of molecule). The left picture shows a phonon 
that shifts the center of gravity of the energy levels, while in the right picture the 
center of gravity is conserved. 

The alkali ions then see only a small change in the net charge and therefore 
couple weakly. 

To analyze the situation more closely, we consider the change in the energy of 
the molecular orbital a under a deformation of the molecule with amplitude 

Q 

Ae° = g, (9) 



where Qa is the electron-phonon coupling. An illustration is given in Figure 2. 
The change of the on-site energy e° will induce a response of the ti„ electrons. 
Since the effect of a point charge q is just to shift the on-site energy hj qU 
(cf. (7)), the screening charge induced by the change Ae° is given by 

Arze. = -7 -jf, (10) 

where 7 > measures the efficiency of the screening: An = —7 q. The total 
screening charge induced by the molecular deformation Q is then given by 

An = ^An„ = -^ {^g^) Q. (11) 

Including screening, the shift in the molecular levels is then given by 

Ae„ = Ae°+C/^An^= L_^^^J g. (12) 
P \ P J 

For molecular solids like the doped FuUerenes the electron-phonon coupling is 
given by the dimensionless constant [14] 
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Given efficient metallic screening (7 1), the coupling to phonons that cause 
a net shift of the molecular levels J2 9/3 ^ will be reduced, while modes that 
leave the center of gravity of the molecular levels unchanged {J2 9/3 = 0) will 
not be affected. Such modes are the Hg modes in Cgo- For these modes efficient 
screening serves to reduce the electron-electron repulsion without affecting the 
electron-phonon coupling. 



6 Conclusion 

We have described how to construct a model for the tiu electrons in the doped 
Fullerenes, that can be analyzed by many body techniques. Using quantum 
Monte Carlo we have calculated the static screening of a point charge. We 
find that the RPA works surprisingly well, almost up to the Mott transition. 
The metallic screening helps to reduce the electron-electron repulsion in the 
doped Fullerenes, where retardation effects are inefficient. But the screening in 
general also reduces the electron-phonon coupling. In a molecular solid there 
can be, however, intra-molecular modes that are not screened, examples being 
the Hg modes in the Fullerenes. We thus find that, although superconductiv- 
ity in the Fullerenes is driven by the electron-phonon coupling, it differs in 
important points from the textbook picture of superconductivity. 
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